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I. INTRODUCTION 
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We study the restoration of chiral symmetry in linear sigma models with two quark flavors. The 
models taken into consideration have a U(2)l X U(2)r and an 0(4) internal symmetry. The physical 
mesons of these models are a, pion, n and ao where the latter two are not present in the 0(4) model. 
Including two-loop contributions through sunset graphs we calculate the temperature behavior of 
the order parameter and the masses for explicit chiral symmetry breaking and in the chiral limit. 
Decay threshold effects introduced by the sunset graphs alter the temperature dependence of the 
condensate and consequently that of the masses as well. This correctly reproduces a second-order 
phase transition for the 0(4) model and for the U(2)l x U(2)r model with an axial U(l) anomaly 
as expected from universality class arguments. Chiral symmetry tends to be restored at higher 
temperatures in the two-loop approximation than in the Hartree-Fock approximation. To model a 
restoration of the axial (7(1) symmetry we imply a temperature-dependent anomaly parameter that 
sharply drops at about 175 MeV. This triggers the restoration of chiral symmetry before the full 
' symmetry is restored and lowers the transition temperatures significantly below 200 MeV. 
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The Lagrangian of masslcss quantum chromodynamics (QCD) with Nf quark flavors has a chiral SU(Nt)i, x 
t"*** ■ SU(Nf)n x U(1)a x U(l)v symmetry. A chiral quark condensate (qq) m (300 MeV) 3 spontaneously breaks the 
\ SU(Nf)A x U(1)a — U(Nf)A part of the symmetry and generates N'j Goldstone bosons. Apart from that there is 
• also a violation of the U(1)a symmetry by instantons 01013 gi ym g ma ss to one of the Goldstone bosons. The U{l)v 
I/-) (vector) symmetry represents baryon number conservation, is always fulfilled and therefore will not be considered 
■ here. Adding mass terms like m q qq to the QCD Lagrangian breaks the symmetry explicitly and gives all Goldstone 
O^l ' bosons a mass, making them pscudo- Goldstone bosons. 

At high temperatures chiral symmetry is expected to be restored. For finite quark masses this happens in a crossover 
transition such that the symmetry is (almost) restored when the temperature (to the third power) is of the order of 
the condensate (qq) 1 / 3 « 300 MeV. Recently, lattice QCD has been able to determine the critical temperature of the 
chiral phase transition. For three flavors it has been found to be in the vicinity of 155 MeV while for two flavors it is 
about 170 MeV for a vanishing quark chemical potential In spite of being a challenging first-principle approach 

to QCD lattice calculation suffer from technical difficulties for small quark masses or for a chemical potential of 
the order of the temperature or larger Q • 

A different nonperturbative approach to QCD is the construction of low-energy effective theories of hadrons with 
' the same chiral U(Nf) x U(Nf) symmetry. The color degrees of freedom are integrated out so that the low-energy 
behavior of QCD is governed by the lightest hadrons which are scalar and pseudoscalar mesons with, in general, light 
quark content. These particles can be found in linear sigma models Q- Since these models have the same symmetry 
as the underlying fundamental theory and thus belong to the same universality class they can be used to study the 
dynamics of phase transitions at finite temperature. Pisarski and Wilczek Q found that there can be a second-order 
phase transition in presence of an explicitly broken U(1)a symmetry, whereas without this axial U(l) anomaly the 
transition is of first order. The two-loop approximation investigated within this article will correctly reproduce these 
features in both the U(2) L x U(2) R and the 0(4) model. 

Linear sigma models cannot be solved analytically so one has to make use of approximations. One problem arising 
at finite temperature is the breakdown of perturbation theory; at a temperature T, a (perturbative) expansion in 
powers of a coupling g yields a new mass scale gT that occurs in the denominators of loop graphs and cancels powers 
of the coupling constant in the perturbation expansion So, terms of all orders of the coupling must 

be taken into account via resummation to avoid these unwanted cancellations. The resummation scheme we a pply 
here is the so-called two-particle point-irreducible (2PPI) effective action introduced by Vcrschclde and Coppcns (13j. 
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Up to the level of the Hartree-Fock approximation it is identical to the two-particle irreducible (2PI) effective action 
formalism by Cornwall, Jackiw and Tomboulis |l4j |. 



Linear sigma models with a U(Nf) x U(Nf) symmetry and two to four quark flavors have been studied in the 
artree-Fock or Hartree approximation within the last more than 25 years [H E [H GJ. Th e O(N) model has 
received even greater attention; it has been analyzed u sing different resummation techniques, where various authors 



received even greater attention: it nas been analyzed using dinerent resummation techniques, wnere various authors 
used local resummations [fj. \wl I2H. I21L I22I I23I I24I I2M l26j |. while, now aday s, nonlocal schemes, like the two-particle 
irreducible effective action |14|. have become popular as w ell l^l^l2^l2ftl3(Tll3lLl32T |. Furthermore, renormalization 




has become a heavily studied issue in this context |33j, |3j> |35j, |36|, |3J 

The Hartree-Fock approximation as well as the two-loop approximation of the 2PPI effective action violate Gold- 
stone's theorem because the formalism's variational parameter associated to the Goldstone boson mass achieves finite 
values at temperatures above zero pfl 123. l2fil |. This problem can be overcome either by looking at the external 
(or physical) propagators — the derivatives of the one-particle irreducible (1PI) effective action — or by a 

construction described by Ivanov et al. |33. For a renormalization group invariant approach to the O(N) model see 
recent work of Destri and Sartirana [3fil4fll|. 

The U(Nf)L x U(Nf)n linear sigma model contains two U(Nj) isospin multiplets — a scalar and a pseudoscalar 
one — each of which is decomposed into an isosinglet and an (Nj — l)-dimcnsional isospin multiplct. For two flavors 
and unbroken isospin symmetry (m u = m^) we obtain four different mesons in the model, a [called Jo (600) nowadays] 
with an isotriplet of (identical) ao bosons in the scalar sector, and r\ with three pions in the pseudoscalar sector. The 
O(N) linear sigma model only consists of a a meson and N — 1 pions and is, for N = N'j = 4, a limiting case of the 
U(2)l x U(2)n model for an infinitely strong U(1)a anomaly. 

This article is organized as follows. In Sees. ITTI and IIIII we describe the U(2)l x U(2)r and the O(N) linear sigma 
model and their pattern of symmetry breaking. Section ll VI deals with parameter fixing and numerical results in both 
models. Finally, in Sec. we draw our conclusions and give an outlook. There is also an Appendix in which more 
details about the computation of the effective action of the U(2)l x U(2)r model are given. 



II. THE 17(2) z, x U(2) R LINEAR SIGMA MODEL 
A. Classical action 

The Lagrangian of the U(Nt)n x U(Nf)z, linear sigma model is given by 

JSf[$] = Tr (d M &&*® - m 2 $ f $) - Ai [Tr ($ f $)] 2 - A 2 Tr [($ f $) 2 ] 
+ c[det $ + dot $ f ] + Tr [H($ + $ f )] . 

The field $ is a complex Nf x Nf matrix containing the scalar and pseudoscalar mesons, 

3? = T a (<j a + in a ). (2) 

Here a a are the scalar fields with J p = + while 7r a denotes the pseudoscalar ones with J p = 0~. The last term in 
the Lagrangian Q describes the interaction with an external field H that breaks the symmetry explicitly, 

H = T a h a . (3) 

T a are the generators of the group U(N f ) such that Tr (T a T b ) = S ab /2. The U(N f ) algebra is fulfilled 

[T a ,T b ] = if abc T c (4a) 
{T a ,T b } = d abc T c (4b) 

where f abc and d abc are the antisymmetric and symmetric structure constants of U(Nf) and a, 6, c = 0, . . . , Nj — 1. 
They are identical to those of SU(Nf) (with all indices starting from one), however for U(Nf) there is in addition 

/ 2 

fabO = , dabO = \ / -Tf 5 ab . (4c) 

V N f 



In the following we will deal only with the case Nf = 2 which reduces the structure constants to 

fijk=£ijk and d ijk = for k € {1, 2, 3} , 



(5) 



where Eijk is the Levi-Civita symbol. The usual identification of the physical bosons for N f = 2 is (see, e.g., Ref. [Tg) 



Since isospin symmetry is left untouched the masses of all particles of one isovector are identical, i.e., m a o = m a ± 
and m-n-o = m^± . 




B. Breaking the symmetry 

The first three terms of the Lagrangian (JIJ are invariant under the group U(2)l x U{2)r ~ 5C/(2)y x S[/(2)ji x 
[/(l)^ x The U(l)v vector symmetry reflects baryon number conservation of QCD. We will not deal with this 

symmetry in this paper as it is always conserved. Chiral symmetry is spontaneously broken if the vacuum expectation 
value of the field $ does not vanish 

(*) = T a 4> a . (7) 

The vacuum should be of even parity, so only c/) a = (a a ) is allowed. According to a theorem by Vafa and Witten 
global vector-like symmetries (isospin, baryon number) cannot be broken spontaneously. So, the remaining symmetry 
must be, at least, SU(2)v- Spontaneously breaking SU(2)a X U(1)a yields four Goldstone bosons, r\ and three pions. 
The determinants in the Lagrangian (JTJ break the U (1)a symmetry explicitly which represents the U (1) axial anomaly 
PJ whose strength is given here by the constant c. This anomaly makes the isosinglet Goldstone boson r\ massive. 
The remaining SU (2) symmetry (of three pions) stays intact if we assume the masses of the up and down quark to 
be equal so that only one diagonal generator gets a finite expectation value. So, we choose 

(*) = T </>o = \cf> 1 . (8) 

Finally, the last term in the Lagrangian explicitly breaks chiral symmetry and makes also the pions massive. It 
resembles the mass terms in the QCD Lagrangian where here H corresponds to the quark mass matrix and $ to the 
quark condensate. We will only deal with the case feo ^ and keep the SU(2) isospin symmetry (m u = mj) conserved 
so that /i 3 = 0. 

With rising temperature we expect the chiral SU(2)v x SU(2)a — SU(2)l x SU(2)r symmetry to be restored so 
that the chiral partners (a and 7T, rj and ao) become degenerate in mass. A violation of the axial U(l) symmetry is 
inherent to the linear sigma model since its strength is directly given by the model's parameter c. The restoration of 
this symmetry can only be modelled in a phenomenological way by making c temperature-dependent, e.g., go down 
with rising T. For c-*0 we expect the rj mass to become identical to the pion mass above a certain temperature so 
that there is a full U(2)a symmetry in the pseudoscalar sector. And, finally, all four masses are expected to become 
degenerate at temperatures above (qq) 1 ^ 3 ~ 300 MeV for explicit symmetry breaking or above a critical temperature 
T c in the chiral limit. 



C. Effective action 

We compute the effective action using the 2PPI formalism 0, 153, H3{ and include all graphs up to two loops. 
The reader is referred to the Appendix for details of the computation. Here, we will only give the final result for the 
effective potential: 

V cS {M 2 , 4> ) = Vd.(<fo) + V dh (M 2 , <f>o) + V q {M 2 , O ) ■ (9) 

It is a function of the masses M 2 , M 2 , M 2 and M 2 o (denoted by M 2 for brevity) and the condensate </>q. The classical 
part is 

VoM = \{™ 2 - c)d? 2 + (^- + ^4- ho 0o . (10) 

The quantum part V q contains all two-particle point- irreducible |6ll | (2PPI) graphs that can be made of the vertices 
of the shifted Lagrangian except for the double bubbles which will be taken care of by Vdb- The propagators within 




FIG. 1: The sunset graphs as in Eq. il l 2t . Scalars are drawn in blue, pseudoscalars in red. Dashed lines represent isotriplets 
while solid lines are used for isosinglets. 



these graphs arc defined by an effective mass and have the Euclidean form 

Here, we only take into account 2PPI graphs with one and two loops which leads to 

V q {M 2 l( f> ) = i lndct(d 2 + Ml) + l - lndet(9 2 + M 2 ) 

+ | lndet(9 2 + M 2 J + | lndet(9 2 + M 2 ) 
+ V SVLnscts (M 2 , (t) ) 
The sunset contribution is given by (note the minus sign) 

V sunmts (M 2 ,(t>o) = -<j>l 



(11) 



Ai + -A 2 ) (3 



+ Ai + 



(12) 



where Sjjfc denotes a sunset graph made of the propagators of the particles i, j and k. A graphical representation 
of these contributions can be found in Fig. ^ These graphs arise from the possible decays of a — > 7T7t, it — > 7777 and 
a — > ?/7r. 

Vdb is the double-bubble part that receives a special treatment in the 2PPI formalism 



y db (M 2 ,0 o ) 
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(3A 2 + 15A 2 o +6A,A ao ) 
Xi 3. 



A CT A„ - 3 



(13) 



A 2 A„„A n 
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A 2 A ff A„ - 3 -Ai + -A 2 A ao A w 



All quantities A arc obtained via 



A«(M 2 ,fo) = 2^V q (M 2 l( f> Q ) (14) 



where * stands for a, ag, r\ or 7r and M 2 is a short-hand notation for all four masses. So, these quantities are explicit 
function of the mass matrix M 2 and the condensate </>o- 

D. Equations of motion 

The mass gap equations in the 2PPI effective action formalism are given by 

dV eS 



= with * = £7, do, T], IT 



and read explicitly 



M 2 a =m 2 - c + 3 ^Ai + y)<& + 3 ( A i + ( A - + A -) 
+ 3 (A! + ^A 2 ) A ao + (a x + A, 
=m 2 _ c + ^ + $ + Ui + ^ (A ff + 5 A ff ) 
+ ^Ai + ^A 2 ^) A„ + ^3 Ax + 7 -\ 2 \ A aa 
M 2 =m 2 + c + (x, + <& + (\x + ^ 3 (A„ + A ao ) 
Ai + y) A ff + 3 (Ai + ^A 2 ) A t 
Ml =m 2 + c + f Ai + ^A 2 ) 4>l + Ui + ^ (A„ + 5 A ao ) 
Vi + |a 3 ) A ct + (^3A 1 + ^A 2 



(15a) 



(15b) 



(15c) 



(15d) 



The formalism is made such that these equations resemble those of the Hartree-Fock approximation with the decisive 
difference that here A* is not a single bubble but calculated from Eq. I|14|) • Neglecting the sunsets contributions to V q 
in Eq. (jl4|) would reduce all quantum corrections to simple bubbles; in this way the Hartree-Fock approximation is 
regained. Solving Eqs. (|T3|) for all four quantum corrections one finds a (quite lengthy) expression for each quantum 
correction in terms of all four masses and the condensate. The equation of motion for the condensate 

^-y eff (M 2 ,0 o ) = o 



can be put in a very easy form by simplifying Vdb{M 2 , <fio) + V C 



ho = Ml O - (2Ai + A 2 ) (g + -£-V sunscts . (16) 



Neglecting the sunsets this equation is equivalent to the one found in the Hartree-Fock approximation [Tlj . 



III. THE O(N) LINEAR SIGMA MODEL 
A. Langrangian 



The linear sigma model with an 0(4) symmetry is obtained from the U(2)y x U(2)a model in the limit of infinite 

2 

0(4)- 



anomaly c — > oo with fixed (m 2 — c) — > ttIq (4 s. The masses of both the 77 and the ag mesons become infinite and thus 



these two mesons drop out of the spectrum. The coupling is renamed to A = (Ax + 4r), and a and the three pions 
now share one 0(4) multiplct. Extending the isospin symmetry from four to TV dimensions we can now write down 
the well-known Lagrangian of the O(N) linear sigma model 

= -\™ 2 <5>} -^($ 2 ) 2 + ^$, with i = l,...,N (17) 

which has been studied extensively 0, EM E3, [2^, H3, E3| ■ In this article we will extend the analysis 

performed earlier j2£j, to the case of explicit symmetry breaking and realistic values of the parameters. 

B. Equations of motion 

The equations of motion are obtained from the 2PPI effective potential where the vacuum expectation value is set 
to be 

= 4>oS Ql , 

so that the 2PPI effective action reads |20L l26l 



Feff^o; Ml Ml) = \ul<& - ±<fi - h 0o + 2A( ^ +2) [Ml + (N- 1)M*] 



1 



(N + 1)M* + 3(N-1)M* (18) 
+ V q (<p ;Ml,Ml) . 



8A (N + 2) 

2 (N - 1) Ml Ml + 2N\ 2 m 4 



The quantum corrections consist of the following one- and two-loop terms 

1. . ,„o ..n, N-l 



V q (<h; Ml Ml) = - lndet(d 2 + M 2 ) + — ^— lndct(d 2 + M 2 ) 



-\ 2 4>l [3 S a<7<7 + {N - 1) S a7r 
The mass gap equations follow from the stationarity conditions 

-V cS = and — — V e g = 



dM a Ln 9M W 
They read 

Ml = m 2 + 3A</» 2 + A[3A (T + (iV-l)A 7r ] (20a) 
Ml = m 2 + A0q + A [A CT + (N + 1) A w ] , (20b) 

where all quantum corrections A are explicit functions of both the condensate 4>q and the masses defined as 

dV 

A i (M a ,M n ;cj> ) = 2-^ ? . 

The equation for the condensate has the same structure as the one in the U(2)l x U{2)r model, Eq. (|16() . 

h = {Ml-2X^l-2\ 2 [?,S <J<J<J + (N -l)S^}} 0o . (21) 
In the following we will only investigate the case N = 4. 



IV. NUMERICAL RESULTS 
A. Parameter fixing and loop graphs at finite temperature 



The parameters in both models are fixed such that at T = the values of all masses are equal to the values in the 
Particle Physics Booklet , cf. Table [I] where we choose the mass of the a meson to be 600 MeV. The value of the 



condensate fa is related to the mesons decay constants f a and determined by the PCAC (partial conservation of axial 
vector current) hypothesis 

fa = d aa Q 00 = 4>0 ■ 

This fixes the condensate to </>o = f n because d aa o = 1, so all decay constants are the same. 

For the two models of this paper the fixing can be done in a unique way because there are as many equations 
of motion as parameters. At tree-level it could be done even in the chiral limit (with ho fixed to zero) since the 
equation for the condensate coincides with the one for the pion mass. Problems occur if one wants to include 
terms that contain a renormalization scale because, first, in the chiral limit there is only one possible value for this 
scale where the parameters can be fixed and, second, for explicit symmetry breaking the temperature-dependence 
of the condensate and the masses is varying with the renormalization scale |Tflj . So, the system only gets an extra 
parameter and all quantities are logarithmically dependent on this scale which makes the results somewhat arbitrary. 
Furthermore, terms orig inatin g from renormalization can be such that the gap equations are not solvable above a 
certain temperature |23l liil Eq . Lcnaghan and Rischke have also shown that, in the O(N) model, there is no 
qualitative difference whether one includes the finite renormalization terms or not. In order to get rid of this extra 
parameter we take the phenomenological approach proposed before [TfiL ITfiL Il9j | and set all finite terms arising from 
regularization equal to zero which makes all quantum corrections only play a role at finite temperature. So, effectively 
the parameters are fixed at tree-level. The resulting values can be found in Tables ITT1 and Hill 

Neglecting finite terms from renormalization the one-loop graphs at finite temperature — the boson determinants 
in the effective action and the single bubble (or tadpole) B — are given by the following equations |2lJ 



lndet(<9 2 +M?) = -5 / dp p 2 In 



1 _ e -E( P )/T 



8 i r°° d 2 



(22a) 
(22b) 



where Ei(p) = \Jp 2 + M 2 and rn is the Bose-Einstein distribution 

f 



n t (p) 



e Ei(p)/T _ I 



(23) 



The sunset graph with three different masses Mi, Mj and Mk is composed of three parts in each of which two of the 
three particles arc taken from the heat bath 



Siik — S, 



»J* — °(ij)k + °(ki) 



Sf'+S, 



(2) 
(jk)i 



(24) 



Here, is a sunset graph with i and j being thermal lines |2£ 
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o «***i^"i(Pi)n*Cft) ^ 



Y + 



Y, 



(ij)k 



(25) 



with 



Y ± = 

(ij)k 



{E. l + E 3 f-{E^ )k f ■ (E t -E 3 f-(E± )k ) 



\ij)k 



and 



E 



(ij)k 



M 2 k 



B. 0(4) model 

For a given temperature T we let N = 4, fix the value of 4>o and then numerically extremize the effective potential 
in Eq. I|18(l with respect to M a and M„. Thereby, we obtain a 1PI potential that is only a function of 4>q. Using the 
equation of motion for the condensate l|21J) we eventually find the temperature-dependent value of the order parameter 
U(T). This quantity is plotted in Fig. El For comparison we also plot the result for the Hartree-Fock approximation 



mass in MeV 
explicit symmetry breaking chiral limit 

particle with anomaly without anomaly with anomaly without anomaly 



a 


6UU.U 


DUU.U 


iiAA A 

6UU.U 


£?AA A 
DUU.U 


7T 


139.6 


139.6 








ao 


984.7 


984.7 


984.7 


984.7 


V 


547.3 


139.6 


547.3 






TABLE I: Meson masses in the linear sigma models with and without axial anomaly for explicit symmetry breaking and in the 
chiral limit. Note that, of course, there is no ao and r\ meson in the 0(4) model. 



explicit symmetry breaking chiral limit 

parameter with anomaly without anomaly with anomaly without anomaly 



m 2 


-(103.78 MeV) 2 


-(388.34 MeV) 2 


-(173.87) MeV) 2 


-(424.26 MeV) 2 


c 


(374.22 MeV) 2 





(387.00 MeV) 2 





/'!) 


(121.60 MeV) 3 


(121.60 MeV) 3 








Ai 


-19.30 


-35.70 


-19.14 


-37.63 


A 2 


78.49 


111.29 


82.73 


119.71 


4>o 


92.4 MeV 


92.4 MeV 


90 MeV 


90 MeV 



TABLE II: Parameters in the U(2)l X U(2)r model for masses as in Table □ 



and for a simplified two-loop approximation (HF in 2-loop) which consists in solving the gap equations (|20(l in the 
Hartrcc-Fock approximation but the condensate equation 1)21(1 with the sunset contributions. In the chiral limit the 
transition temperature is almost the same as in the true two-loop approximation though the condensate drops faster 
at lower temperatures and then eventually approaches zero at T w 210 MeV. In the chiral limit of the Hartree-Fock 
approximation, there is a first-order phase transition at T c ks 181 MeV whereas the two-loop approximations (both 
the true and the simplified one) correctly reproduce a second-order transition though at a higher temperature of about 
210 MeV. 

For explicit symmetry breaking the simplified two- loop approximation results lie on top of the Hartree-Fock results, 
only the true two-loop approximation yields slight deviations. There, the condensate drops faster at temperatures 
below 200 MeV and decreases more slowly for temperatures above 250 MeV than but in all approximations the 
crossover temperature — the one where the slope of f-w(T) is largest — is about 225 MeV. 

The temperature dependence of the a and pion mass is displayed in Fig. [3] In the chiral limit [Fig. |3(a)| the a mass 
behaves similarly in both approximations while the pion mass first increases before it slightly drops at T w 195 MeV 
(180 MeV in the chiral limit). Beyond temperatures of about 200 MeV the pion mass in both approximations is almost 
the same. In the vicinity of 300 MeV both masses become identical, a sign for the restoration of chiral symmetry. 
There, the temperature (to the third power) is equal to the value of the chiral quark condensate (see above). In the 
chiral limit [Fig. |3(b)| the pion mass first grows stronger in the two-loop approximation than in Hartree-Fock, then 
it slightly drops before the critical temperature is reached and the masses become the same. 

Moreover, there is a violation of Goldstone's theorem at finite temperature in the chiral limit because even for 
temperatures lower than the critical one, M n is not equal to zero although the symmetry is spontaneously broken. 
This phenomenon can also be found in ear lier results on the O(N) model d . We will comment on this 
in Sec. 



parameter 


explicit symmetry breaking 


chiral limit 


m 2 


-(388.34 MeV) 2 


-(424.26 MeV) 2 


ho 


(121.6 MeV) 3 





X 


19.94 


22.22 




92.4 MeV 


90 MeV 



TABLE III: Parameters in the 0(4) model for masses as in TableQ 
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FIG. 2: Temperature dependence of the condensate in the 0(4) model. Comparison of the two-loop approximation with the 
Hartree-Fock approximation (HF). 



C. U{2) L x U(2) R model 

The procedure performed for the O(N) model, cf. Sec. lIVBl cannot be done for the U{2)l x U(2)r model because 
it turns out that the potential V = V c \ + Vdb has only a saddle point with respect to the four masses instead of a local 
extremum. 

Solving the gap equations (|15f) is a cumbersome procedure because they contain derivatives of the sunset graph 
with respect to a mass. And in the vicinity of the decay threshold of the particles involved in a sunset graph, e.g. 
M a ~ 2 M ff , the sign of the derivative of the sunset graph quickly changes (see Fig. which results in a numerically 
unstable behavior in this region. To avoid this trouble a simplified approximation is used (cf. Sec. IIVB|l : we solve 
the mass gap equations <|15|) in the Hartree-Fock approximation, i.e., the quantum corrections A to the masses 
only consist of single bubbles given by Eq. (|22b(l . We substitute these masses into the two-loop equation for the 
condensate Ijl6(l to obtain the temperature-dependent order parameter or decay constant Jtt{T) and, for simplicity, 
call this approximation "two-loop" from now on. 



1. Model with axial U(l) anomaly 

In Fig.[S]we show the temperature dependence of the condensate in the U (2)l x U(2)r model with an axial anomaly. 
In the chiral limit of the Hartrcc Fock approximation there is a first-order phase transition at T c ps 178 MeV where 
the condensate discontinuously drops down from about 52 MeV to zero. In the two-loop approximation we find 
a second-order transition at T c as 272 MeV. So, the two-loop approximation correctly reproduces the order of the 
phase transition obtained from universality class arguments (3|, though the critical temperature is about 100 MeV 
higher than suggested by QCD lattice calculations 0, 0] . The artificial first-order transition in the Hartree-Fock 
approximation has been discovered before by Roder et al. for a a mass of 400 MeV . For both explicit symmetry 
breaking and in the chiral limit the order parameter decreases more slowly with rising temperature in the two-loop 
approximation than in the Hartree-Fock approximation and exhibits some "bumps" in the curve. The reason for 
that behavior is a decay threshold effect caused by the sunset graphs [cf. Eq. 1)12(1 and Fig.0]. To check where the 
thresholds are crossed we plot the mass ratios for the decays a — ► 7T7t, ao — > Tpr and a — > rji] in Fig. EJ 
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(a) Explicit symmetry breaking. 
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FIG. 3: Temperature dependence of the masses in the 0(4) model. Comparison of the two-loop approximation with Hartree- 
Fock (dotted). 
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FIG. 4: The ratio S n7T7T /S aaa for M a = 500 MeV at T = 150 MeV and T = 300 MeV, and M a = 400 MeV at T = 200 MeV 
and T = 250 MeV. 



The curve of the order parameter in the two-loop approximation deviates from the one of the Hartree-Fock ap- 
proximation at a temperature of about 175-180 MeV. This is the region where the thresholds of the decays a — > tttt 
and ao — + rjir are crossed (see Fig. EJ. The sunsets become larger with rising temperature but also decrease with 
masses approaching the threshold from below. The latter behavior suddenly changes at the threshold and the sunsets 
suddenly grow which causes the immediate deviation from the Hartree-Fock results. Looking at the equation for 
the condensate <|16f) one can roughl y co nclude that a larger value of the condensate is needed to compensate for the 
contribution from the sunset terms |62j . 

Comparing the prefactors of the different sunsets in Eq. H12J1 using the numerical values of the couplings from 
Table ITT1 we notice that, in the case of a finite U(1)a anomaly, the contributions from S avv and S aoV7T are almost of 
the same size whereas S amr is smaller by a factor of eight. Furthermore, the sunset term has dimension two and thus 
scales with temperature and the masses of the particles within. So, considering the prefactors and the different values 
for the masses — where M ao and M v are the two largest ones — the deviation from the Hartree-Fock approximation 
is dominated by a threshold effect of the decay ao — > ryK. 

The temperature-dependent masses for the U(2)l X U(2)r model are displayed in Fig. Especially the masses 
of the scalar mesons a and ao behave differently in the two-loop approximation than in Hartree-Fock whereas the 
masses of both pseudoscalar mesons exhibit no qualitative difference in their temperature dependence between the 
two approximations. The observed deviations are due to the aforementioned threshold effects through sunset graphs. 
At a temperature of about 300 MeV (or beyond T c w 272 MeV in the chiral limit) the masses of the chiral partners 
become identical so that chiral symmetry is restored. The U(1)a symmetry remains broken since the parameter c 
is not a function of temperature. The gap between the mass squares of the isospin partners r\ and 7r (or a and ao) 
remains equal to 2c according to Eq. i|15[l . In the U(2)l x U(2)r model as well, the pion mass is finite in the chiral 
limit even for T <T C which was discovered in the Hartree-Fock approximation earlier . We will give a statement 
concerning a possible violation of Goldstone's theorem in Sec. 



2. Model without axial U(l) anomaly 



For the U(2)l x U{2)r model without axial U(l) anomaly we plot the condensate vs. temperature in Fig. El In 
the chiral limit there is a first-order phase transition at T c ps 200 MeV (T c w 170 MeV in Hartree-Fock) which is 
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FIG. 5: Temperature dependence of the condensate in the U(2)l x U(2)r model with an axial U(l) anomaly. Comparison of 
the two-loop approximation with Hartree-Fock (HF) for explicit symmetry breaking (ESB) and in the chiral limit. 



expected from the respective universality class [jj. 

Here, the deviation from the Hartree-Fock approximation sets in at lower temperatures because the a — > 7777 
threshold is crossed already at temperatures about 100 MeV, followed by a — > mr at about 175 MeV and, finally, 
the on-shell decay ao — > rjir becomes impossible at 200 MeV (see Fig. EJ. Compared to the case with a finite axial 
anomaly the effect of the sunset Sa^* is only about 7% (instead of 12%) of that of the other two. So, in this case as 
well, the deviation from the Hartree-Fock approximation is dominated by S ao riTT- 



In Fig. 10(a) we observe that, for explicit symmetry breaking, the masses of the chiral partners tend to become 
identical at 250 MeV (200 MeV in Hartree-Fock) bef ore the actual restoration of the full U(2)l x U(2)r symmetry 
takes place at about 300 MeV. In the chiral limit, [Fig. 10(b)| these two points coincide when the critical temperature 
is reached. At finite temperature the mass of the 77 meson differs from the pion mass although they both started from 
139.6 MeV (or zero in the chiral limit) at zero temperature. This indicates that the approximations considered in this 
article are not well-suited to model the 77 meson as a fourth (pseudo-)Goldstone boson since they seem to contain an 
effective ^(1)^ breaking through the une qual treatment of 77 and 7r in the gap equations l|TH]l. This has been found 
before in the Hartree-Fock approximation (l5|. In the case of a vanishing U(l) anomaly as well, Goldstone's theorem 
seems to be violated in the chiral limit as mentioned before because M v and M n arc finite below T c . 



3. Temperature- dependent U(1)a anomaly 

There are indications from the lattice that at high temperatures effects arising from the U(1)a breaking are strongly 
suppressed 0, 0> 0, |5(1 HH1 EE EE 0- This suggests an effective restoration of the U(1)a symmetry close 
to the critical temperature. We model this by fixing the parameters at zero temperature to the physical masses but 
considering the anomaly parameter c as a function of temperature. As an example we describe a suddenly dropping 
behavior at 175 MeV with the Fermi function 

C ( T ) = T- rr, WAT11 (26) 

1 + cxp[(T - T A )/AT\ 

with T A = 175 MeV and AT = 10 MeV, see Fig. El The values are chosen such that at T sa 200 MeV the strength 
of the anomaly has dropped by one order of magnitude. 
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FIG. 6: Thresholds for the decays a —* tttv, ao — * t/tt and a — > r\r\ in the U(2)l X U(2)r model with axial anomaly as. 
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FIG. 7: Masses in the U(2)l x U(2)r model with axial anomaly as functions of temperature. Comparison of the two- loop 
approximation with Hartree-Fock (dotted lines). 



> 

CD 



100 
90 
80 ■ 
70 ■ 
60 ■ 
50 ■ 
40 ■ 
30 ■ 
20 ■ 
10 ■ 
■ 



2-loop chiral limit 
HF chiral limit 
2-loop ESB 
HF ESB 



50 100 150 200 250 300 350 400 

T/MeV 

FIG. 8: Temperature dependence of the condensate in the U(2)l X U(2)r model without an axial anomaly. Comparison of the 
two-loop approximation with Hartree-Fock (HF) for explicit symmetry breaking (ESB) and in the chiral limit. 



This causes a sharp decrease of the order parameter at about 175 MeV, cf. Fig.^] In the chiral limit, there is a very 
weak first-order transition at T c = 180.49 MeV where the condensate drops to zero from a value of 1.44 MeV. In the 
Hartree-Fock approximation we find a first-order transition at T ss 161 MeV. Shifting Ta influences the order of the 
transition in the two-loop approximation; for Ta < 175 MeV the phase transition is of first order because the anomaly 
has effectively vanished at T c , and for Ta > 175 MeV there is a second-order transition because the axial anomaly is 
sufficiently strong at the critical temperature. Changing AT in Eq. 1)26(1 has the same effect. For explicit symmetry 
breaking, there is a crossover at a temperature of about 180 MeV. The masses of the chiral partners become identical 
at about 200 MeV (see Fig. ll3fl . In the two- loop approximation, these temperatures are significantly lower than those 
for the models with a fixed U(1)a anomaly. Although the anomaly parameter tends to zero the full symmetry is 
only finally restored at about 300 MeV where all four masses become identical. The temperature behavior of the 
condensate and the masses are quite similar to those found in the framework of a Nambu-Jona-Lasinio (NJL) model 
with three quark flavors [5^, [5(| . 

The effect that a suddenly dropping anomaly parameter triggers the restoration of chiral symmetry was observed 
earlier in the linear sigma model with three quark fla vors Il7l . In the chiral limit, the difference between the scalar 
and the pseudoscalar meson mass beyond T c [cf. Fig. 13(b)] is only driven by the function c(T), Eq. So, chiral 
symmetry restoration is triggered by the effective restoration of the U(1)a symmetry in the chiral limit, too. 



V. CONCLUSIONS AND OUTLOOK 



We have investigated the restoration of chiral symmetry in the U(2)l x U(2)r and in the 0{N) linear sigma models 
at finite temperature. Taking into account two-loop sunset contributions makes the condensate drop less rapidly 
with increasing temperature; the crossover and critical temperatures arc significantly higher than in the Hartree-Fock 
approximation which is shown in Table llVl The deviation from the Hartree-Fock approximation is induced by decay 
threshold effects, especially of the decay ao — > tjtt, which drive the value of the condensate away from zero. For 
the U(2)l x U(2)r model with an axial U(l) anomaly the masses of the chiral partners become identical at about 
300 MeV for explicit symmetry breaking or at a critical temperature of about 270 MeV in the chiral limit, whereas a 
mass gap of 2c between the isospin partners remains. But even with a zero anomaly parameter, there is an effective 
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FIG. 9: Thresholds for the decays a — » 7T7t, ao — > r/n and cr — > 7777 in the U(2)l X U(2)r model without axial anomaly. 
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FIG. 10: Masses in the U(2)l x U(2)r model without axial anomaly as functions of temperature. Comparison of the two-loop 
approximation with Hartree-Fock (dotted lines). 
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FIG. 11: Modelled temperature dependence of the anomaly parameter c. The behavior is like the function in Eq. 1261 with 
T A = 175 MeV and AT = 10 MeV. 
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FIG. 12: Temperature dependence of the condensate in the U(2)l x U(2)r model without an axial anomaly. Comparison of 
the two-loop approximation with Hartree-Fock (HF) for explicit symmetry breaking (ESB) and in the chiral limit. 
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FIG. 13: Masses in the U(2)l X U(2)r model with a temperature-dependent axial anomaly as in Fig. 1111 Comparison of the 
two- loop approximation with Hartree-Fock (dotted lines). 



model 




transition temperature (order of phase transition) 




Hartree-Fock 


two-loop 


0(4) 




FSR 99^ MpV 


99 "i MpV 




rhiral limit 1 K1 MpV 

L.111I cLL illlllLr lOl Ivlc V I lab 1 


91 fl MpV 




with anomaly 


ESB 210 MeV 


270 MeV 


U(2) L x U{2) R 


chiral limit 178 MeV (1st) 


272 MeV (2nd) 


without anomaly 


FSR IQ 1 ! MpV 
jjjO-Lj iyj ivic v 

chiral limit 170 MeV (1st) 


94-^ MpV 
200 MeV (1st) 




with varying anomaly 


ESB 180 MeV 
chiral limit 161 MeV (1st) 


185 MeV 
181 MeV (weak 1st) 



TABLE IV: Transition temperatures and order of the phase transitions of the two models in two different approximations. For 
explicit symmetry breaking there is always a crossover. 



U(1)a breaking by the approximation itself due to an unequal treatment of the 77 meson and the pions in the gap 
equations l|15f) . Here as well, the sunsets only raise the temperature where chiral symmetry is restored compared to 
the Hartree-Fock approximation. 

We have also investigated the effect of a temperature-dependent anomaly parameter c(T) as in Eq. I|2f>[l. With a 
steep decrease at about 175 MeV, such that the strength is reduced to 10% at temperatures of about 200 MeV, we 
could reduce the transition temperature to 181 MeV in the chiral limit of the two-loop approximation. The phase 
transition is (very weakly) of first order for this parameters, though it can be changed to second order by letting c 
decrease slower with rising temperature. The masses of the chiral partners become identical at significantly lower 
temperatures (about 200 MeV for explicit symmetry breaking) than with a fixed anomaly parameter. Nevertheless, 
the full symmetry is also only restored at about 300 MeV (explicit symmetry breaking) or at 220 MeV where both 
the condensate and the anomaly have vanished. Though it is obvious that a suddenly dropping anomaly parameter 
triggers the restoration of chiral symmetry in the two-loop approximation. 

There se ems to be a viola tion of Goldstone's theore m in the chiral limit because all masses of the Goldstone bosons 
[M n in Fig. |7(b) and 13 (b)| and M v in Fig. 10(b)| are finite below the critical temperature T c . The reason for that 



is the fact that resummation schemes usually violate global symmetries |3J|. The physical masses, obtained as second 
derivatives of a one-particle irreducible effective potential, are not identical to the variational mass parameters. The 
usual geometric argument that proves the physical Goldstone masses to be zero ist the following. Consider the second 
derivative of a 1PI effective potential depending on the O(N) invariant cf) 2 = (\> a <i>a 



d 2 



= 2V'{$ 2 )5 ab 

={$) 



4V"{P) </>*<!>„ 



<P=(4>) 



(27) 

$=($) 



With the expectation value pointing only in the 0-direction, (cf>) = (0o, 0, . . . , 0), the Goldstone masses are given by 
the derivatives perpendicular to that direction and thus are proportional to V'{4> 2 ) which is zero for a non-trivial 
vacuum. So, whenever there is a minimum off zero in the potential there are Goldstone bosons |63j . 

Comparing this work with recent publications [23, [2~3| one has to state that the effect of non-local corrections to 
the propagators seems to be more efficient for the reduction of the critical temperature from the Hartree-Fock value 
of about 200 MeV to a desired value of about 175 MeV 0, |f| than considering only local corrections. The approach 
described in this article is can only reduce the crossover temperature significantly if a varying anomaly strength is 
taken into account. So, it would be interesting to apply non-local resummation schemes also to the U(2)l X U(2)ji 
model. 

Furthermore, the effective U(1)a symmetry breaking that is inherent to the Hartree-Fock and two-loop approxi- 
mation may be remedied by a different approximation, possibly by one inspired by a 1/iV expansion similar to that 
used in Ref. j^. 

Including strange mesons (Nj = 3) could possibly lead to interesting non-linear effects since the U(1)a anomaly 
term is trilincar for three flavors and thus would generate additional sunset graphs with different signs. Adding 
fermions (nucleons or constituent quarks) to the model using a Yukawa coupling is also an attractive extension of this 
work |57j . For the O(N) model this has been done by several authors before 
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APPENDIX A: EFFECTIVE ACTION OF THE U{2) L x U(2) R MODEL 



This Appendix contains details for the computation of the effective action of the U(2)l x U(2)r model. We define 
the vacuum expectation value 



- a *Va 



(3) = T a 

with a real-valued condensate <f>o [cf. Eq. (jHJ] and shift the (complex) U(2)l x U(2)r fields to 

$= ^(p 1 + T a (a a + iTT a ) , 

where a a and 7r a are real and symbolize the scalar and pseudoscalar meson fields. The shifted Lagrangian is a sum of 
four parts, 

J2f[*] - -Ki as8 (0o) + + + ^[*] ■ (Al) 
The first part is the classical tree-level potential 

Vd ass (0) = 2 m2< ^o _ 3 ^ ab ^ a< ^ b + g-^afccd <t>a<i>b§c<t>d ~ h a (j) a , (A2) 

which for the choice of Eq. (jHJ reduces to the expression in Eq. I)10[l. The structure of the mass matrices and 
interactions is given by the coefficients |l5j 

Sab = - (^a0<5o0 - ^al^fcl - <5 Q 2<5b2 - ^3(^3) (A3a) 

Fabcd = "J" (^ab<5cd + <^ad^6c + <5ae<5bd) 
A 2 

dnbc + dacndnbd) (A3b) 

o 

Habcd = ~ p^ab^cd H J~ (d a bnd n cd + facnfnbd + fbcnfnad) ■ (A3c) 

The second part of the shifted Lagrangian consists of all bilinear terms in the fields 



%l\$a.\ CT a , 7T„] = TjidfiCTa) 2 - - [m 2 S ab -QQab + ^Fabcd </>c</>d] ^"a ^b 



(A4) 



+ 2 ~('^ 7r a) - 2 I-" 1 <5ab + 6 ^ afc + 4 ^ abcd TTa^b 



Looking at the structure of this equation we conclude that the tree-level mass matrix is diagonal if the expectation 
value <f> a has only one component, e.g., <p a = 5 a o<po- So there is neither a mass mixing among the scalar and 
pseudoscalar mesons nor between them. According to the 2PPI resummation scheme 0, |2(J, [2^, [2(j one constructs 
the one-loop terms as given in Eq. . 

The third part of Eq. IjAlJ) contains trilinear terms describing the three-particle vertex 

4 

■&3[<f>ai °~a, 7T a ] = - ^ J^abcd ^d^a^b^c ~ ^Habcd 4>d^a^b^c ■ ( A 5) 

The sunset contributions to the effective potential arc constructed from that by computing the expectation value 

V sunsct (M 2 ,<f )0 ) = -l( [ [ J? 3 (x)J? s (y)\ (A6) 



2 



2PPI 



where we chose conventions for Euclidean space-time and neglected an overall volume factor on the left hand side. 
The subscript 2PPI means that only two-particle point-irreducible graphs may be constructed from this expectation 
value resulting in sunset graphs with three propagators from the vertex at x to the vertex at y as required by Wick's 
theorem. Note that since there is no tree-level mixing of masses, valid contractions are only of the type 

(<r a (x) a a {y)) — > 1 or (tt (x) ir a (y)) — > 9 1 (A7) 

p z + M* a p 2 + M 2 a 

where <j$ = a, 01,2,3 = &o an d fto = 77, ""1,2,3 = tt- This finally yields the sunset contribution to the effective potential 

+ 2^ 2 ^ ^ Scree + 3 S aao ao + S<nm) 

2 3 1 (A8) 

+ ( 2" j Sa-rrir + 2^ 2 ^OofT 

where Syfc denotes a sunset graph with propagators of the particles i, j and fc. The last part of the shifted Lagrangian 
contains the rest, i.e., the four-vertex interactions, 

-24[c ,7r ] = - 2H a bcd OaOb^c^d - T^Fabcd {o~ a (J b (T C U d + ^a^b^c^d) ■ (A9) 

We can now write down the general structure of the two-particle point-irreducible (2PPI) effective potential 

V e g(M 2 , 0o) = F class (0 o ) + 2 W Qbc , Af b Af d + ^ Qbc(i (Af b Af d + Af 6 Af d ) + V q (M 2 ,<j> ) , (A10) 

where 

&db = {<ra(x)<T b (x)) and Af 6 = (7r (x)7r 6 (x)) (All) 

denote expectation values of local composite operators ( "bubbles" ) , and V q contains higher-order corrections made of 
2PPI graphs likes sunsets, basketballs etc. For the masses are diagonal the expectation values (|A11(I are diagonal as 
well. Using the physical identification of the meson fields as given in Eq. 10 we declare 

A a = A^ A„ = AS, (Al2a) 

A O0 = A s n = Af a = Af 3 A w = Af x = A^ 2 = Af 3 . (A12b) 

Contracting the coefficients T a b c d and Tiabcd with the appropriate expectation values we obtain the following expression 
for the 2PPI effective action 

V cS {M 2 , 0o ) - Klass^o) + V q (ct>o, M 2 ) 
Ai A 2 



, (3 A 2 + 15A^ + 6A CT A 7r ) 

Y + Y ) ( 3 A * + 15 < + 6 A v A -o ) (A13) 



' 4 ; ■ - v 2 '4 

y + ^A 2 ) A, A, + 3 0Ai + ^A 2 ) A O0 A„ , 

where T4g is a function of the condensate (j>o and all four masses M a , M ao , M v and M„. Note that this is the case for 
each quantity A* as well since it is given by Eq. (|14|l . Equation (|A13|) looks like the double-bubble term in the 2PI 
expansion [l5j except for the sign which is different here due to the special role of this contribution and the quantity 
A within the 2PPI formalism pi I2I l2i | . 

For a detailed graphical comparison of the 2PPI and the 2PI effective action see the Appendix of Ref. [23] . 
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